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When a laser pulse impinges on a molecule which is invariant under certain symmetry opera- 
tions selection rules for harmonic generation (HG) arise. In other words: symmetry controls which 
channels are open for the deposition and emission of laser energy — with the possible application 
of filtering or amplification. We review the derivation of HG selection rules and study numerically 
the interaction of laser pulses with an effectively one-dimensional ring-shaped model molecule. The 
harmonic yields obtained from that model and their dependence on laser frequency and intensity are 
discussed. In a real experiment obvious candidates for such molecules are benzene, other aromatic 
compounds, or even nanotubes. 



I. INTRODUCTION 

Harmonic generation (HG) from atoms (L'Huillicr & Balcou 1993), molecules (Liang et al. 1994), clusters (Donclly 
ct al. 1996), and solids (von der Linde et al. 1995) as a short-wavelength source is of great practical relevance. In 
recent years a huge amount of publications were devoted to "harmonic engineering," under which we subsume cither 
the study of phase matching during the propagation of the emitted light through gaseous media (Gaarde et al. 1998), 
generation of attosecond pulses using varying ellipticity (Antoine et al. 1997), multi-color studies (Milosevic et al. 
2000), or HG from thin crystals (Faisal & Kaminski 1996 and 1997). Experimental results on HG from cyclic organic 
molecules were also reported in the literature (Hay et al. 2000a and 2000b). However, the authors made no attempt 
to verify theoretically predicted selection rules. In the well established physical picture of HG one assumes that an 
electron tunnels out of the atom or ion, moves in the laser field and eventually rescattcrs with its parent (or other) ion 
where it might recombine — leading to an emission of a photon with several times the fundamental frequency (Becker 
ct al. 1997). Assuming such a viewpoint one can easily explain prominent features of HG spectra such as the famous 
cut-off at I p + 3.17U P in the single atom case (where I p is the ionization energy and U p is the ponderomotive energy, 
i.e., the cycle-averaged quiver energy of the electron in the laser field). For a general overview on harmonic generation 
in laser fields see the recent review by Salieres et al. (1999). 

Fortunately, to derive the selection rules for harmonic emission only symmetry considerations are necessary (Alon 
ct al. 1998). We now briefly summarize this approach. Let us consider a Hamiltonian which is periodic in time, i.e., 
H(i) = H(i + t). Such a Hamiltonian might describe an electron in a long laser pulse where the pulse envelope is 
sufficiently adiabatic. The time-dependent Schrddingcr equation (TDSE) reads (we use atomic units K = e = m = 1 
throughout, if not noted otherwise) 



* £ (r,t)=0. (1) 



According to the Floquet theorem (see, e.g., Faisal 1987) we can write ^>£(r,t) = ips(r, t) exp(— if t) with ip£(r,t) = 
ipsir^t + t) leading to the Schrddinger equation 

H F {t)^ s (r,t) = £Mr,t), H F {t) = H{t) - i- (2) 

at 

with £ the so-called quasi energy and Hf the Floquet Hamiltonian. Since ^£(r,t) is periodic in time it might be 
expanded in a Fourier series ipsir, t) = J2 n ex P[~ in(u>t + S)]^ 1 ^ (r). This Floquet approach is a well-known method 
in multiphoton physics (Gavrila 1992) and a widely used numerical simulation technique also (Potvlicgc 1998). In 
order to derive the HG selection rules we now assume (for reasons which will become clear soon) that the system is 
in a single and non-degenerate Floquet state ip£. The HG spectra (HGS) peak no. n is present only if the Fourier 
transformed dipole moment pt(r) does not vanish for the frequency nu), 
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J dt exp(-ino>i) J d 3 r^* £ (r,t)fj l (r)^ £ (r,t)^0. (3) 

Now let us suppose we worked out a symmetry operation V under which the Floquet Hamiltonian is invariant, 
VHf^V^ 1 = H_f(<). It follows that (in the case of non-degenerated Floquet states ips, see assumption above) 
Vipe = dips holds. Here a is a phase factor, \a\ = 1. For convenience, we rewrite ([|) as 

((* £ |)Lt(r)exp(-inw*)|*£)) = ((ip £ \p(r) exp(-inut)\tpe)) ^ (4) 

where the double brackets indicate spatial and temporal integration [this is the so-called extended Hilbcrt space 
formalism, see, e.g., Sambe (1973)]. The HG selection rule can be derived from 

({ft>e\n(r)exp{-irujt)\il/e)) = {{V^ £ \T^{r) eKv{-inut)V- l \Vip £ )) 

= ((#|pM(r)exp(-ina;t)p- 1 |#)) (5) 

leading to 

/i(r) exp(— inuit) = Vp(r) exp(— \nutyp~ 1 . (6) 

The last step in (||) is not possible if ip £ is not a pure Floquet state. In the next Section we apply (^|) to derive the 
selection rule for HG from ring-shaped molecules. Here, for the sake of illustration, we rederive the selection rule 
for a single atom with spherically symmetric potential V(r) and linearly (in x-direction) polarized laser field E(t) = 
Ee x smut. The Floquet Hamiltonian Hp(t) = — ^V 2 + V(r) + Exsmut — idt is invariant under the transformation 
x — ► —x and t — ► t + n/u. If we look for harmonics polarized in the x-direction we have fJ>(r) = x and from (|^) 
xexp(— inut) = — xexp[— inu(t + tt/u)] follows exp(— imr) = — 1. Therefore only odd harmonics are generated. 

The paper is organized as follows: in Section |h| we introduce our model and the selection rule which holds in its 



case. In Section [II we present and discuss our numerical results. Finally, in Section IV we give a summary and an 
outlook. 

II. A SIMPLE ONE-DIMENSIONAL MODEL FOR RING-LIKE MOLECULES 

The time-dependent Schrodinger equation (TDSE) for a single electron in a laser field E(t) and under the influence 
of an ionic potential V{r) reads in dipole approximation and length gauge 



i-¥(r, t) = ^--V 2 + V(r) + E(t) ■ rj *(r, t). (7) 

The dipole approximation is excellent since in all the cases studied in this paper the wavelength of the laser light is 
much greater than the size of the molecule. 

If we force the electron to move along a ring of radius p in the zy-plane V(r) becomes V(tp) where <p is the 

usual polar angle. With an electric field of the form E(t) = E(t) £,cos(ut)e x + \fl~- ~ C 2 " sm{ut)e y , where E(t) is a 

slowly varying envelope, u is the laser frequency, and £ is the ellipticity parameter, the TDSE (0) becomes effectively 
one-dimensional (ID) in space and reads 



q fid 2 r l\ 

i— V((p,t) = i-^-2g-^ + V{tp)+E(t)p [£ cos(ut) cos ip + y/1 - e sm{ut) sin J t) . (8) 

In case of circularly polarized light (£ = 1/ v2) this simplifies to 

i^*(v,*) = f~^+nv) + ^oofl( V -«t)J*( Vl t). o) 

We now assume that the potential V(ip) has an TV-fold rotational symmetry, V(ip + 2ir/N) = V(<p). Then, with the 
help of (Q), we can easily derive the selection rule for HG in the system described by the TDSE (||). The transformation 
(tp — > tp + 2%/N, t — > t + 2n/Nu) leaves the corresponding Floquet Hamiltonian invariant. For, e.g., anti-clockwise 
polarized emission p,(r) = pexp(ii^) holds, and from (||) we have 
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pexp(ic/?) exp(— inujt) — pcxp[i(ip + 27r/iV)] cxp[— inu;(i + 2tv /Nu>)] 

leading to n = Nk + 1, k = 1,2, 3.... For the clockwise emission one finds accordingly n = Nk — 1. Thus we expect 
pairs of HG peaks at kN ± 1 (Alon et al. 1998). 

The TDSE (|) or (§) can be easily solved ab initio on a PC. We did this by propagating the wavefunction in 
time with a Crank-Nicholson approximant to the propagator U(t + At,t) = exp[— iA£H(i + At/2)] where H(t) is the 
explicitly time-dependent Hamiltonian corresponding to the TDSE (g). Our algorithm is fourth order in the grid 
spacing Atp and second order in the time step At. The boundary condition is ^(O,*) = ^(2^, t) for all times t. 

III. NUMERICAL RESULTS AND DISCUSSION 

We now present results from single active electron (SAE)-runs with p = 2.64 (bond length and radius of benzene 
CeH 6 ) and an effective model potential 

^) = -y[cos(7V^) + l] (10) 

with N = 6 and Vq — 0.6405. This leads to an electronic ground state energy £o = —0.34 which is the experimental 
ionization potential for removing the first electron in benzene (see, e.g., Talebpour et al. 1998 and 2000). Note, that 
in our simple model we have no continuum but discrete states only. The first six excited states are located at —0.27 
(two-fold degenerated), —0.07 (two-fold degenerated), 0.16 (non-degenerated), 0.48 (non-degenerated), 0.85 (two-fold 
degenerated), 1.48 (two-fold degenerated). The energy levels of our model resemble, apart from an overall downshift 
and the removal of degeneracies of certain states, those of the isoperimetric model where Vq = 0, the energy levels 
are given by £ m = m 2 /2p 2 , m = 0, 1, 2, . . . with the states m ^ two-fold degenerated. Therefore, the energy level 
spacing, and thus typical electronic transitions, are different from those in real benzene. However, it is not our goal 
to present quantitatively correct results for laser benzene-interaction in this paper but we rather want to demonstrate 
some of the underlying principles of HG from ring-shaped molecules in general. 

In Fig. [l] we present HGS for a q = 240 cycle pulse of the shape E(t) = Esm 2 (ujt/2q) with an electric field 
amplitude E = 0.5 a.u. and frequency cu = 0.18. In Fig. [l]a the result for linear polarization £ = 1 is shown, in Fig. [If} 
the result for circular polarization £ = l/v2- To obtain those plots we evaluated the Fourier transformed dipole in 
x-direction, i.e., ((^(t)\pcos((p) exp(— inut)^ (t))) . In the £ = 1-case the dipole with respect to y is clearly zero while 
in the circular case there is simply a phase shift of it/2 with respect to the dipole in x. As expected, in the linearly 
polarized field all odd harmonics are emitted whereas in the circular case only the harmonics 6fc ± 1, k = 1,2,... are 
visible. Other emission lines arc many orders of magnitude weaker. Those lines can form band-like structures which 
are interesting in itself. However, in this paper we focus only on the laser harmonics. They dominate the HGS, at 
least as long as no resonances are hit. 

In Fig. H the emitted yield of the fundamental and the first four harmonics (5th, 7th, 11th, 13th) in a circularly 
polarized laser pulse are presented as a function of the laser frequency. The pulse length T = 8378 a.u. (corresponding 
to sw 200 fs) and peak field strength E = 0.2 a.u. (corresponding to 1.4 x 10 15 W/cm 2 ) were held fixed. The frequency 
is plotted in units of the smallest level spacing of the model molecule, i.e., the energy gap between first excited and 
ground state, which is in our case Q = 0.34 — 0.27 = 0.07. For laser frequencies w < fi we do not find the 7th or 
higher harmonics. The 11th and 13th harmonic show an overall decrease with increasing laser frequency whereas 
the fundamental, the 5th and the 7th stay relatively constant in intensity. For frequencies u> < 2.50 there is a 
complicated dependency of the harmonic yield on ui. All harmonics show a local maximum around 1.351. However, 
at that frequency one apparently hits a resonance since the harmonic peaks become broad and show a substructure 
(see left inlay in Fig. ||). In the interval 2 < u)/fl < 2.5 the fundamental drops whereas the 7th harmonic increases 
in strength. Note that the 7th harmonic is anti-clockwise polarized, like the incident laser field, whereas the 5th is 
polarized in the opposite direction. For frequencies u> > 2.5Q, the behavior becomes more smooth apart from another 
resonance near 3.8f2. In general, the HGS look clean for sufficiently high frequencies and far away from resonances 
(like in the right inlay of Fig. ^). A rich substructure near resonances is visible in the HGS (cf. inlay for to = 2fl). 

In Fig. U harmonic yields for the fixed frequency u = 2.8Q as a function of the field amplitude E are shown. The 
higher harmonics (llth-19th) appear only for higher field strengths whereas fundamental, 5th and 7th are rather 
weakly dependent on the field strength. The anti-clockwise polarized harmonics (polarized like the incident laser 
light, drawn thick) tend to overtake the clockwise polarized ones (drawn thin) at higher laser intensities. However, 
E = 0.6 corresponds already to a laser intensity 1.3 x 10 16 W/cm 2 where a real benzene molecule would probably 
break. 
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It is interesting to study the scaling of the TDSE (||) with respect to the size of the molecule. If one scales the 
molecule radius like p' = ap, the TDSE (||) remains invariant if t' = a 2 t, V' = V/a 2 , E' = Ejo? . u>' = uj/a 2 is chosen. 
From this and our numerical result that laser frequencies > fl are preferable for clean HGS wc learn that molecules 
bigger than benzene are more promising candidates for HG with realistic laser frequencies [from Nd (to = 0.04 a.u.) 
to KrF (cu = 0.18 a.u.)]. 

One might object that electron correlation could spoil the selection rule because in reality it is not only a single 
electron which participates in the dynamics. However, Alon et al. (1998) have proven that this is not the case. This is 
due to the fact that (i) the electron interaction part of the Hamiltonian is still invariant under the transformation V, 
and (ii) V commutes with the (anti-) symmetrization operator. Even approximate theories or numerical techniques 
like time-dependent Hartree-Fock or density functional theory do not spoil the selection rule since they involve only 
functionals which depend on scalar products of single particle orbitals, all invariant under the transformation V 
(Ceccherini & Bauer 2001). 

IV. SUMMARY AND OUTLOOK 

In this paper we demonstrated numerically HG from a model molecule with discrete rotational symmetry which is 
subject to a circularly polarized laser pulse. In particular the harmonic emission from an effectively ID model with 
-/V = 6 (i.e., a simple model for benzene) as a function of laser frequency and intensity was discussed. It was found 
that for frequencies below the characteristic level spacing fi HG is strongly affected by resonances. The situation 
relaxes for higher frequencies. For the efficient generation of higher harmonics laser frequencies > Q and rather strong 
fields are necessary. In such fields real aromatic compounds probably ionize and dissociate already. 

Numerical studies for a more realistic, effectively 2D model molecule will be presented elsewhere (Ceccherini & 
Bauer 2001). 

In order to obtain short wavelength radiation it is desirable to have either k or N in the selection rule kN ± 1 as big 
as possible. From our numerical results we infer that it is probably hard to push efficiently towards high k without 
destroying the target molecule. For that reason, in a real experiment nanotubes are promising candidates because N 
can be of the order of 100 or more (Dresselhaus et al. 1998), and, moreover, HG should be even more efficient when 
the laser propagates through the tube. However, it remains the problem of the proper alignment of the laser beam 
and the symmetry axis of the molecule. Crystals might be better candidates in that respect. 
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FIG. 1. HGS for a 240 cycle sin 2 -shaped pulse with an electric field amplitude E = 0.5 a.u. and frequency ui = 0.18. The 
polarization was (a) linear (£ = 1) and (b) circular (£ = 1/^/2). In (a) odd harmonics up to n — 23 dominate whereas in (b) 
the harmonics obey the selection rule 6fc ± 1, k = 1, 2, . . ., i.e., the 5th, 7th, 11th, 13th, 17th, 19th harmonics are visible. 

FIG. 2. Fundamental and harmonic yield vs. laser frequency lo in units of the smallest level spacing Q = 0.07. Laser 
field intensity and pulse length was kept fixed. The anti-clockwise polarized fundamental, 7th, and 13th are plotted thick, 
the clockwise polarized 5th and 11th are drawn thin. The inlays show spectra (harmonic yield vs. harmonic order) for three 
particular frequencies (indicated by arrows). See text for further discussion. 



FIG. 3. Harmonic and fundamental yields for fixed frequency lo — 2.8Q but different field amplitude E. The inlays show 
spectra (harmonic yield vs. harmonic order) for three particular E (indicated by arrows). See text for further discussion. 
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Fig. 1: D. Bauer and F. Ceccherini, "A numerical ab initio study ..." 




Fig. 2: D. Bauer and F. Ceccherini, "A numerical ab initio study ..." 




0.20 0.30 0.40 0.50 0.60 

E [a.u.] 



Fig. 3: D. Bauer and F. Ceccherini, "A numerical ab initio study ..." 



